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Numerical study of liquid metal flow in a rectangular duct 
under the influence of a heterogenous magnetic field 

Evgeny V. Votyakov^ Egbert A. Zienicke 



Abstract: We simulated numerically the laminar flow 
in the geometry and t he magnetic field of the experi- 



menta l channel used in | Andreev. Kolesnikov. and ThessI 
(|2006r )]. This provides detailed information about the 
electric potential distribution for the laminar regime (nu- 
merical simulation) and in the turbulent regime as well 
(experiment). As follows from comparison of simulated 
and experimental results, the flow under the magnet is de- 
termined by the interaction parameter N = Ha^ /Re rep- 
resenting the ratio between magnetic force, determined 
by the Hartmann number Ha, and inertial force, de- 
termined by the Reynolds number Re. We compared 
two variants: (i) (/?f',A^)=(2000,18.6) (experiment), 
(400,20.25) (simulation), and (ii) {Re,N) =(4000,9.3) 
(experiment), (400,9) (simulation) and found an excel- 
lent agreement for the numerical and experimental dis- 
tributions of the electric potential. This is true despite 
of the fact that the experimental inflow is turbulent while 
that in the simulation is laminar. As a special feature of 
the electric potential distribution local extrema under the 
magnets are observed, as well experimentally as numeri- 
cally. They are shown to vanish, if the interaction param- 
eter falls below a critical value. Another interesting new 
detail found in our numerical calculations is the appear- 
ance of helical paths of the electric current. Using a sim- 
plified magnetic field without span-wise dependence, we 
show that important physical features of the considered 
problem are sensitive to small variations in the spatial 
structure of the magnetic field: the local extrema of the 
electric potential and also the helical current paths dis- 
appear when the simplified magnetic field is used. The 
structure of the three dimensional velocity field is also 
investigated, in particular, a swirling flow is found in the 
comers of the duct caused by Hartmann layer destruction 
behind the magnets. 
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1 Introduction 

The flow of an electrically conducting fluid under a lo- 
calized inhomogeneous magnetic field is of interest for 
many industrial applications deahng with the problem 
to influence hot metal melts by the use of magnetic 
fields. This has the advantage that no direct co ntact with 
the ch emically aggressive melt is necessary I Davidsor] 
(Il999l) l. One prominent example is the electromagnetic 
brake us ed in modem continuous steel casting, see for in- 
stance iTakeuchi. Kubota. Miki. Qkuda. and Shiroyama 
(l2003l) 1. If one neglects the liquid steel jets entering 
the mould, the flow in a liquid metal channel under a 
static localized magnetic field may reproduce qualita- 
tively many features of the electromagnetic brake: the 
braking effect on stream-wise velocity, suppression of 
turbulence under the magnet, and effects of strong spa- 
tial dependence of the magnetic field. Another exam- 
ple, where liquid metal flow in a channel is important for 
possible industrial applications, is the Lorentz Force Ve- 
locimetry based on exposing the fluid to a magnetic field 
and measu ring the drag force acting upon the mag netic 
field lines liThess. Votyakov. and Kolesnikovl (120061) 1. 

In the recently appeared experi mental work 
I Andreev. Kolesnikov. and Thess the liquid 



metal flow was systematically investigated for the range 
of interaction parameters A <N <20. In the experiments 
special attention was focused on the suppression of 
turbulence by the magnetic field and on a systematic 
recording of data for the electric potential building up 
under and around the magnets. As a result, the full map 
of the electric potential distribution was obtained for the 
middle horizontal plane of the rectangular duct and few 
A'^ parameters. If the transverse magnetic field would be 
homogenous the electric potential data might be used to 
determine the veloc ity components perp endicular to the 
magnetic field LSommeria and Moreau (fl 982) 1. In the 
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Figure 1: Coordinate system and sketch of the prob- 
lem: rectangular channel 2Lx x 2Ly x 2Lj and two mag- 
nets on the bottom and top walls. Qualitatively, there 
are shown transverse magnetic field B^{x) varying along 
stream-wise direction, and span-wise profile of stream- 
wise velocity Ux{y). The center of the coordinate system 
is in the center of the magnetic gap. 



present case of strongly inhomogeneous magnetic field, 
the electric potential data can be solely compared with 
the results of the corresponding numerical simulations. 

Fig- [H presents schematically the geometry 
and magnetic field config ur ation studied in 
l Andreev. Kolesnikov. and ThessI (2006)]: the liq- 
uid metal moves in a rectangular duct, where the locally 
heterogenous magnetic field is created by two permanent 
magnets on the top and bottom walls of the duct. The 
originally convex velocity profile Uy{x) adopts, by 
passing the magnetic field B^{x), a characteristic M 
shape what is a manifestation of the electromagnetic 
brake process. 

The main goal behind our numerical investigation is first 
to reproduce features of the flow and the electromagnetic 
quantities found in the experiment. Moreover, since nu- 
merical calculations have the advantage that all quantities 
are available in the whole computational domain, the sec- 
ond goal is to visualize the additional data, i.e. 3D veloc- 
ity field and electric current paths, which are not acces- 
sible to measurements in the experiment. Before we ex- 
plain why one can successfully compute the experimen- 
tal flow under and near the magnet using laminar calcu- 



lations let us shortly describe in the next paragraphs the 
state of the art of numerical and theoretical approaches 
that could be appropriate to describe the present problem 
involving turbulence and inhomogeneous magnetic field 
at the same time. 

The most ideal approach from the point of view of phys- 
ical exactness would be to carry out a direct numerical 
simulation (DNS) of the experimental flow which mod- 
els correctly the whole ifow including turbulent regions 
down to small scales. However, to do this under a mag- 
netic field steeply varying with space coordinates is ex- 
tremely hard since it requires space and time resolution 
a few orders of magnitude more than is available today. 
Actually, to catch all spatial structures of the flow, espe- 
cially Hartmann layers and sidewall jets, the simulation 
has to be fully three-dimensional and needs high reso- 
lution near the boundaries. This concerns also the tur- 
bulent flow which was generated in the experiment by 
a honeycomb at the inflow in order to observe suppres- 
sion of fluctuations by the static magnetic field. The in- 
flow distance takes almost half the length of the compu- 
tational domain, therefore half of the computational re- 
sources must be paid to simulate correctly the flow with- 
out explicit action of the magnetic field, just to monitor 
declining turbulent fluctuations. 

Other numerical approaches to catch turbulent features 
of the flow are Large eddy simulations (LES) and 
Reynolds averaged stress models (RANS). These are 
as well faced with serious technical difficulties. The 
present state of the art for LES of conducting fluids 
under magnetic fields is such that one may treat the 
flow inside a homogenous magnetic field only if one 
resolves the boundary lay ers as in a direct numeri - 
cal simulation (DNS), see | Knaeoen and Moin (|2004h : 



Vorobev. Zikanov. Davidson, and K naepenI ( 2005 )1. The 
definition of appropriate wall functions for the turbu- 
lent Hartmann and sidewall layers is still an actual 
field of research. LES in inhomogeneous magnetic 
field meets serious difficulties as one has to find ap- 
propriate methods of spatial averaging. For the mo- 
mentary available RANS models it is not yet clear 
whether they are able to describe the unavoidable 
anisotropy of the tur bulent scales inside strongly varyin g 
magne tic field, see jWid lund. Z ahrai. and BarkI (Il998h : 
Kenjere_candHanjalic| (2000, 200^7 In any case, before 
to start a LES or RANS study one first has to define and 
verify parameters appearing in these phenomenological 
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models by DNS runs. 

The usual analytical means also fail to describe the fea- 
tures of the experimental flow. This holds, because the 
parameters of the applied magnetic field lie outside the 
limits imposed by the assumptions necessary for any an- 
alytical theoretical treatment. In particular, to neglect in- 
ertial effects, the existing analytical approaches assume 
very strong and slowly varying magnetic field, which is 
not fulfilled neither for the electromagnetic brake nor 
for the channel flow considered here. Typically, any 
regular theoretical method is based on an asymptotic 
expansio n of MHD equations around large [see for 



examp le iLavrentiev. Molokov. Sidorenkov. and Shishko 
(Il990h and references therein]. However, in the sys- 



tem under consideration, the local interaction number 
A'^(x) = B^ {x)N goes up from zero to moderately high val- 
ues on a short distance under the inward gradient of the 
magnetic field. Moreover, even if such an approach is not 
entirely impossible for some cases, it employs a series 
expansion what is of the same cost approxi mately as a 
full 3D simulation iSellers and Walkeil (ll999h l. Also, the 
theoretical methods do not take a span-wise dependence 
of the magnetic field into account, however, as we shall 
see later, this seemingly fine detail of the magnetic field 
configuration can be responsible for significant qualita- 
tive features in the electric potential distribution inside 
the magnetic gap. 

However, as follows from the exp erimental data 
I Andreev. Kolesnikov. and ThessI ( 2006 )1. the intensity 
of turbulent fluctuations inside the magnetic gap is lower 
than one percent, and essentially smaller than at the inlet 
distance at front of the magnetic system. This provides 
evidence that the magnetic field is strong enough to be 
the main influence shaping the flow structure inside the 
magnetic gap. This flow structure, qualitatively charac- 
terized by a M-shaped profile, is weakly dependent on the 
separate Re and Ha numbers and strongly dependent on 
the interaction parameter N, especially when A'^ is high. 
Another conjecture from these experimental data is that 
the originally turbulent inlet velocity profile is of small 
importance as well. 

The foregoing statements give us an opportunity to re- 
produce the experimental results using a laminar numer- 
ical 3D simulation. That means, we do not carry out a 
computation with the same values of Re and Ha as in the 
experiment and do not monitor turbulence. Instead, since 
the main effects are due to the interaction number A'^, one 



may select for the simulation Re and Ha lower than in the 
experiment but belonging to the same ratio N = Ha^ /Re. 
As inlet flow one takes a laminar duct flow. To clarify, 
whether the shape of inflow velocity plays a role for the 
measured experimental data, one can test different lam- 
inar inflow profiles having the same mean flow rate but 
different flatness. 

It is easy to point out a reasonable range for the Re 
numbers in the numerical simulation. Large Re param- 
eter provokes turbulence which could not be properly re- 
solved with current computational resources while too 
low Re number results in a viscous force in the core of the 
flow. On the other hand, the highest limit of the Re num- 
ber is governed also by the value of Ha number which is, 
in its own turn, dictated by the available grid resolution 
for the Hartmann layer. 

Thus, the main goals of the simulation were to determine 
the qualitative global velocity field and to find a good re- 
production of the electric potential in the magnetic field 
region. Especially, we were interested to find the two 
extrema of the electric potential which were observed in 
the experiment. As it turns out, all these aims were reach- 
able. The overall general features of the laminar flow are 
well represented by our code. In section 14.11 we present 
how the Hartmann layer and the sidewall jets are form- 
ing under the magnet in a stationary flow. For the elec- 
tric potenti^ljinde£_the^^ experimental results 
of IIAndreev. Kolesnikov. and ThessI (12006) 1 show clearly 
(section I4i2l) that the flow under the magnet is near to be- 
ing laminar. In this region the magnetic and the inertial 
forces are predominant, if one excludes the regions very 
near to the walls where the viscous forces are essential. 
Therefore we find a good representation of the electric 
potential distribution of the experiment for our runs us- 
ing lower Hartmann and Reynolds numbers but keeping 
the same interaction parameter. 

The structure of the present paper is the following. In 
section [2] the equations and our numerical method to 
solve them are presented. Here also the inflow profiles 
and the used grid are specified. In section |3] some proper- 
ties of the experimental magnetic field are explained and 
a second simplified magnetic field with no span-wise de- 
pendence is introduced, which serves to show that small 
changes on the inhomogeneous magnetic field can lead to 
remarkable differences in the electric potential distribu- 
tion. In section HI we present the results of our numerical 
computations showing all characteristics of the velocity 
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field (section |4~TI ). tlie comparison of the numerically de- 
termined electric potential with experimental data (sec- 
tion |4]2l), and the electric current paths (section |43] ). 

2 Equations and numerical method 

The governing equations for electrically conducting and 
incompressible fluid are derived from the Navier-Stokes 
equation coupled with the Maxwell equations for mov- 
ing medium, and also using the Ohm's law. We apply 
the quasi-static (induction-less) approximation where it 
is assumed that an induced magnetic field is infinitely 
smal l in cornparison to the external magnetic field (see, 
e.g. iRobertsI (Il967h l). so it is neglected when one calcu- 
lates the Lorentz force, but it is not neglected at finding 
the electric current density j. The resulting equations in 
dimensionless form are then given as follows: 



3u , 

3^ + (u-V)u 
V u 
j 

A(|) 



-Vp + ^Au + NiixB), (1) 
Ke 

0, (2) 

-V(|) + uxB, (3) 

V-(uxB). (4) 



Here u denotes velocity field, B is an external mag- 
netic field, j is electric current density, p is pressure, (|) 
is electric potential. Re = uqH /v is Reynolds number, 
= Ha^ /Re is the interaction parameter (Stuart num- 
ber), and Ha = HBo{a/fj)^^^ is Hartmann number, all 
defined with the half-height of the channel H, mean ve- 
locity uq, typical magnetic field strength Bq, density p, 
electric conductivity a, kinematic v and dynamic = pv 
viscosities. 

As follows from eq. ([T]) the viscous force Au is scaled by 
Reynolds number Re, therefore at high Re and far from 
the walls it plays a minor role. As a result the flow is 
governed by the interaction parameter N defining the ra- 
tio between magnetic and inertial forces. This is a case 
we treat in the paper under consideration by comparing 
experimental and simulated results with similar N and 
different (and high) Re numbers. 

In the experiments | Andreev. Kolesnikov. and ThessI 

(12006) 1 ■ the eutectic alloy Gao.68l1io.20Sno.12 was used as 
a hquid metal. It has density p = 6360 kg/m^, elec- 
tric conductivity a = 3.46 • 10^ Ohm^' and kinematic 
viscosity v = 3.4 • 10^^ m^/s. Thus, the Hartmann 
number Ha = HBo{o/pvy^^ defined with half-height 
of the channel {H = I cm) and magnetic field intensity 



So = 0.504 T is Ha = 193. The interaction parameter 
= Ha^ /Re was varied in the experiments by means 
of the mean velocity rate uq entering the Reynolds num- 
ber Re = Huq/v. We have implemented the experimen- 
tal range of interaction parameters by varying mainly Ha 
(up to 120) and keeping Re = 400, thus in the simulation 
4<A^<36. 

Unknowns of the equations ([T] - |4l) are the velocity 
vector field u{x,y,z), and two scalar fields: pressure 
p{x,y,z) and electric potential ^{x,y,z)- The domain 
of the flow is given by a rectangular channel (Fig. [T]) 
(l-^l < l-x, \y\ < l-y, and |z| < with Lx = 25, Ly = 
5, and = H = 1) having the same aspect ratios 
(Length: Width:Heig ht=25 :5 : 1 ) as the ex perime ntal chan- 
nel of I Andreev . Kolesnikov. and ThessI ( 2006 )1. (In this 
experimental paper, = 2 is defined as a whole height 
of the channel, and in the present paper we take // = 1 as 
a half-heig ht. Thus, the Re and Ha nu mbers given in 
llAndreev. Kolesnikov. and ThessI (|2006l ')1 were divided 
by factor two here.) For the external magnetic field we 
used the magnetic field which was measured in the ex- 
periment and a second configuration, which is used for 
comparison (see section O. The Hartmann number is 
based on the Bo value in the center of the magnetic gap 
x = 0,y = 0,z = 0. Note, this So value is not the maximal 
one, to fulfill curl- and divergence free requirements, the 
field is slightly increasing by approaching top and bottom 
walls at fixed x = and y = 0. 

To solve the partial differential equations ([T] - H]) ini- 
tial and boundary conditions have to be provided. Since 
we are interested in a stationary solution the initial con- 
ditions play no role (except for the speed of conver- 
gence), and for the boundary conditions we suppose a 
duct with electrically insulating and "no-slip" walls on 
the sides, top and bottom. Insulating and "no-slip" con- 
ditions require at the boundary Y to impose that M|r = 0, 
3(|)/3n|r = 0, where n is normal vector to T. The outlet 
was treated as a force free (straight-out) border for the 
velocity. The electric potential at inlet and outlet borders 
was taken equal zero because the inlet and outlet are suf- 
ficiently far from the region of magnetic field. 

For the inlet profile we made use of two possibilities. A 
self-consistent choice for laminar flow is the stationary 
laminar profile of an infinite rectangular duct known an- 
alytically in the form of a series expansion. In the ex- 
perimental channel, however, the inflow is generated by 
a honeycomb shaping a more flat inflow profile in span- 
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Figure 2: XY and YZ projections of the grid for the parti- 
tion A^;, X A^j, X A^, = 32 X 32 X 32 and /? = 1 .25, S = 0.75, 
T = 0.75. 



32^, R=1 .25 S=0.75 T=0.75 
64^, R=1 .25 S=0.75 T=0.75 
64^, R=1.25 S=0.85 T=0.85 
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Figure 3: span-wise profile of the stream- wise velocity 
computed at different grid parameters: A'',^^ x Ny x = 
32 X 32 X 32 (squares) and 64 x 64 x 64 (crosses and 
dots), {R,S,T) = (1.25,0.75,0.75) (squares and crosses) 
and (1.25,0.85,0.85) (dots). Other parameters are com- 
mon: Re = 100, A^ = 16, ;c = 0, z = 0. 



wise and vertical direction and generating vortex struc- 
tures which give rise to decaying turbulence on the way 
of the liquid metal to the magnet. To study the influ- 
ence of a more flat turbulent like profile on the electric 
potential under the magnet, we have generated a second 
inlet profile in the following well defined way: periodic 
boundary conditions (outflow=inflow) were imposed in 
stream-wise direction and then turbulent runs for the 
Reynolds number of the experimental system were per- 
formed. Then, by a space averaging procedure we com- 
puted the turbulent mean inflow profile. In spite of this 
careful check of the influence of inlet boundary condi- 
tions, we have found no difference in the final results 
between the laminar and the more flat turbulent inflows, 
except for the transitional region before the magnet. This 
shows that the magnetic field is sufficiently strong to 
completely govern the hydrodynamics of the flow when 
the magnetic region is reached. It agrees also with the 
experimental observation about measured turbulent fluc- 
tuations essentially suppressed inside the magnetic gap. 

As a base for our solver we used NaSt3DGP 
the simulation code developed in 
IIGriebel. Domseifer. and Neunhoeffen (il995 ))1. Orig- 



inally, this finite difference solver was designed for 
pure hydrodynamical problems. Therefore we had to 
extend it by the following features to be able to solve 
MHD problems: (1) using the Poisson solver also for the 
determination of the electric potential, and (2) including 



the Lorentz force contribution on the right hand side 
computing the preliminary velocity field. Moreover, we 
reorganized the input and output parts of NaSt3DGP in 
order to work with the arrays keeping magnetic field, 
electric potential, and electric current. 

Briefly, the numerical algorithm is the following. To de- 
couple Eg. dH) - dS), the Chorin-type projec t ion m ethod is 
apphed jPevret and Tavloil (Il983h : iHirschI (Il988h l. This 
is a general procedure based on a predictor-corrector ap- 
proach. First, the Poisson equation for the electric po- 
tential ([U) is solved and the electric current is found ac- 
cording to ([3]l. Next, a preliminary velocity field is com- 
puted from the momentum equation without regarding 
the influence of the pressure term. The second part of 
the time step from the preliminary velocity field, which 
is not divergence-free, to the new divergence-free veloc- 
ity field allows to derive a Poisson equation for the new 
pressure. Thus, the whole algorithm is written down as 
follows (index n denotes time integration step, the spatial 
discretization is omitted): 

1. Solve Poisson equation for the electric potential: 

= V-(u"xB). 

2. Compute electric current: 

f = -Vf' + u"xB. 
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(a) at z=0 



(c) at y=0 
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(b) B^ at y=0 
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(d) Byaty=-5 




Figure 4: Cuts of the experimental magnetic field: transverse (plot a at z = and plot b aty = 0), stream-wise (plot 
c at J = 0) and span-wise (plot d at y = —5) components. For plot a: y = (solid), y = 3 (dashed), and y = 5 
(dot-dashed); for plots b — d: z = — 1 (solid), -0.66 (dashed), and (dot-dashed). Plots b and c corresponds also to 
the simplified spanwise-homogeneous magnetic field discussed in the text. Bold dashed lines are the borders of the 
physical magnet. 



3. With the known u" and j", find the preliminary ve- 
locity field u 

= —Au"+N{r X B) - (u" • V)u". 
ot Re 

4. To compute the velocity field u"+^ of the next inte- 
gration step one has to solve the Poisson equation 
for pressure 

and, as a result one finds 

u"+^ =u-btVp"+\ 

The above algorithm is explicit and for simplicity it is 
presented as a scheme of first order precision. It de- 
scribes the principal sequence, which in the code is re- 
alized with the Adams-Bashforth time step having sec- 
ond order precision. For pressure stabilization we apply 



a staggered grid, and the VONOS (variable-order non- 
oscillatory scheme) scheme is used to discretize the con- 
vective and diffusive terms. 

The solver is implemented to support parallel computa- 
tion: the channel in the program code is subdivided into 
domains, and each domain is run as a separate process. 
Communication between the processes takes place at ev- 
ery integration step to synchronize the borders between 
the domains. 

To make sure that all the layers are properly resolved in 
the simulation, we use an inhomogeneous 3D grid con- 
structed in the following way. First, we map —L^ <x< 
Lx, —Ly < y < Ly, —Ly < z < L^, onto three auxiliary 
variables — 1 < r < 1,-1 <s<\,—\<t<\a?, follows: 

tKRD tg{S^) tg{T^) 

' th{R) ' ' tg{Sli/2y ' tg{Tii/2y 

and then variables r,s,t are uniformly partitioned into 
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Nx, Ny, parts. There are three stretching parameters 
R,S,T providing a denser (R) grid at x = and close 
to walls {S,T). Typical values used in the simulation 
are R = 1.25, and 0.75 <S < 0.95, 0.75 <T < 0.95. 
Fig. |2] gives an example of the grid for the partition 
Nx X Ny X A^, = 32 X 32 X 32 and R = 1.25, 5 = 0.75 
T = 0.75. For each value of the Hartmann and Reynolds 
number in calculations, parameters R,S, T as well as the 
number of grid points were varied in order to obtain grid- 
independent results, see Fig. [3] The typical values used 
for the simulation were grid 64 x 64 x 64 and R = 1.25, 
5 = 0.85 r = 0.85. 

3 Magnetic Field 

The magnetic field used in the present paper was 
measured at equidistant points in t he experiments 
I Andreev. Kolesnikov. and ThessI ( 20061) 1. including all 
three components, and then interpolated on the grid 
points and stored into a three-dimensional array which 
is supplied as an input to the solver. This is the most 
general approach which enables us to work with any 
magnetic field configuration supplied externally. The de- 
tails about the field are given in t he exp erimental paper 
I Andreev. Kolesnikov. and ThessI ( 200^ 1. here we just 
remind that it is created by means of two permanent mag- 
nets fixed outside on the top and bottom walls of the 
channel as is shown in Fig. [T] In the gap between the 
magnetic poles the field is aligned mainly along the verti- 
cal direction parallel to the z-axis. Outside of the gap, be- 
fore and behind the poles, there are regions characterized 
by an inward and an outward gradient of the transverse 
magnetic field, see Fig. \^a, b), where also the stream- 
wise component, B^, plays a role (Fig. HJc)). Since the 
physical magnets are finite in the span-wise direction (y- 
coordinate), the magnetic field is also dependent on the 
y coordinate, in particular, on the side walls a By com- 
ponent different from zero appears, Fig.|4l;<i). Moreover, 
the vertical component B^ is rather lower near the side 
walls (3^ = ±Lj,) in relation to B^ in the center Cv = 0) 
(Fig. IHa)). A few distinguished cross-sections of trans- 
verse, stream-wise, and span-wise magnetic field compo- 
nents along the x-axis are shown in Fig. HI Other detailed 
plots of the r nagnetic field can be found in the exp eri- 
mental paper I Andreev. Kolesnikov. and Thess (200a)]. 
In addition to the experimental magnetic field configura- 
tion, which changes slightly along the j-axis, a few vari- 
ants were simulated for a simplified magnetic field (in 



the following always referred to as 'simplified magnetic 
field' in contrast to the 'experimental magnetic field') 
being independent of the y coordinate. The span-wise 
homogeneous field in this case originates from external 
magnets which are infinitely long in j-direction. It is 
characterized only by stream- wise Bx{x,z) and transverse 
B^{x,z) components which coincide with Bx{x,y,z)\y=Q 
and B^{x,y,z)\y=Q in the vertical slice {y = 0) shown in 
Fig. Ht^jC). It turns out that the span-wise decline of B^ 
near the side walls as well as that of the other B compo- 
nents is of decisive influence on both the electric poten- 
tial distribution and electric currents under the magnetic 
poles (see discussion in sections |4!2] and 1431 ) . This illus- 
trates an importance of fine details of the localized het- 
erogenous magnetic field, which are usually neglected in 
theoretical computations. 

4 Results of the simulation 

This section is divided in three subsections focusing on 
different physical quantities — the velocity field, the 
electric potential, and the vector field of the current den- 
sity — giving together an overview of the general and 
specific features of the experimental system. The known 
general features are the formation of a Hartmann profile 
in the transverse direction and the formation of an M- 
shaped profile in span-wise direction as well in the re- 
gion of increasing and in the region of decreasing mag- 
netic field (counted in stream-wise direction). In the first 
subsection we demonstrate that these two processes lead 
to a complex three-dimensional flow structure, which in 
the given magnetic field and channel proportions only 
can be determined numerically. For the electric potential 
distribution a direct comparison with experimental data 
is presented in the second subsection. A special feature 
showing up is a pair of extrema of the electric potential 
under the magnet. These extrema disappear when either 
the interaction parameter N is too small or the simplified 
span-wise uniform magnetic field (at any value of A'^) is 
used. The same sensible dependence of phenomena on 
the magnetic field structure is found in the last subsection 
for complicated helical current paths, which are present 
for the experimental magnetic field and are absent for the 
simplified magnetic field. 

4.1 3D velocity field 

We start from the well known phenomenon for duct flow 
under a locally heterogenous magnetic field: the M- 
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x=0 



x=0, z=0 




Figure 5: 3D M-shaped velocity profile, x = 0,Re = 400, 
N = 36 



shaped velocity profile. It is shown in Fig. [5]for Re = 400 
A'^ = 36 at X = 0. Such a profile is a consequence of the 
braking effect of the magnetic field applied to the elec- 
trically conducting and moving fluid. Shortly, the effect 
can be understood in the following way: under the action 
of the externally applied magnetic field, electric currents 
are induced in span-wise direction. The larger the inten- 
sity of the magnetic field B, the higher the density of the 
electric current j. The channel walls are insulating and 
the magnetic field is locally heterogenous, therefore, to 
make a closure the electric currents will leave the area 
of the high magnetic field (see also the figures in sec- 
tion |43]). Then, as j and B are both present, the Lorentz 
force Fi = j X B, hampers stream-wise fluid movement 
in the bulk of the channel. The flow tries to circumvent 
the area of high magnetic field as much as possible, and, 
as a result, the stream-wise velocity profile will adopt M- 
shaped in span-wise direction, and a stagnant region with 
stream- wise velocity about zero (see solid lines in Fig.[6l) 
forms inside the magneticgap. Various M-shaped sur- 
faces are given in iSterll (Il990'l1. including their discus- 
sion and corresponding references. It is worth to note 



that the external magnetic field selected in [Sterl! (119901 )] 
for the simulation was either divergence- or curl-free, but 
not simultaneously divergence- and curl-free, neverthe- 
less this did not disturb the formation of a M-shaped pro- 
file. 

Fig. |6] shows the influence of the interaction parameter 




Figure 6: Stream-wise velocity profiles in span-wise di- 
rection in the center of the magnetic gap (x = and 
z = 0), Re = 400 and = 4 (dot-dashed), = 9 (dashed), 
N = 20.25 (dotted), A'^ = 36, experimental field (solid line 
1) and A'^ = 36, simplified field (solid line 2). 



A'^ on the M-shaped velocity profile. The effect is clear: 
the higher the value of A'^ the stronger the Lorentz force 
braking the flow, the lower the velocity in the center, and 
thus the more pronounced the stagnant region. The solid 
curves 1 and 2 in Fig.[6]are related to the same parameter 
values Re = 400, A'^ = 36, but to different magnetic fields: 
the experimental field (curve 1) and the simplified field 
(curve 2). This comparison shows that the character of 
the flow in this two cases is rather different. In particular, 
for the case of simplified field, the width of the stagnant 
region is larger but its level is higher than in the case 
of the experimental field. This can be explained by the 
decrease of in span-wise direction for the experimen- 
tal magnetic field. As the transverse magnetic field near 
the side walls consequently for the experimental field is 
lower than in the center, the liquid can more easily flow 
around the "magnetic obstacle" (see below), therefore the 
redistribution of the flow for the experimental magnetic 
field is more pronounced than for the simplified field. 

Fig. ID 

Another insight for the M-shaped profile formation might 
be taken from the XY vector velocity plot shown in 
It was noticed several times a kind of similarity be- 
tween the flow of an electrically conducting liquid pass- 
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Figure 7: Vector plot of the velocity field in the central 
horizontal plane, z = 0. The upper part shows the inten- 
sity of the magnetic field on the straight line z = 0,y = 0, 
Re = 400, = 36. 



ing through an area of high local heterogenous magnetic 
field (magnetic obstacle ) and the well known flow around 



a blu ff body (see e.g. jCuevas. Smolentsev. and Abdou 



^06)]). One can see that the velocity vectors envelop 
the central part of the channel nearly in the same way 
as it would be for a solid cylindrical obstacle in a two- 
dimensional flow. The only difference is that a real body 
creates a region with no fluid, say velocity zero, while in 
the region with magnetic field there still is fluid, but with 
very small velocity. 

To analyze span-wise and vertical velocity components 
we have plotted a series of vector plots in vertical sec- 
tions of the channel at several fixed values of x. They 
are shown in Fig. [8] to highlight the intrinsic three- 
dimensional nature of the velocity field. In the region 
of increasing magnetic field before the magnet gap (see 
Fig. the velocity vector in the vertical section is 
nearly aligned towards the comers of the channel. This 
can be explained by the simultaneous action of two ef- 
fects: (i) the beginning tendency to form a M-shaped 
profile, and (ii) the formation of Hartmann layers which 
is accompanied by a flattening of the velocity profile in 
vertical direction (see also Fig.|9]below). The second re- 
organization process of the flow is accomplished faster, 
i.e. finished already at smaller x-coordinate, than the 
first one. This can be concluded from the second sec- 



tion. Fig. [8^7, because the plotted vectors indicate only a 
movement in the span-wise direction. 

The formation of M-shaped profile continues further 
when the flow passes the maximum of the magnetic field, 
see Fig. [8]:. However, as the flow reaches the region of 
decreasing magnetic field, there is an inversion of the 
span- wise movement (Fig. [8]^): now it is opposite to the 
one observed before the magnet. These two processes 
— flow towards the side walls at the front of the mag- 
netic gap and towards the center (i.e., away from the 
side walls) after the gap — demonstrate again that the 
flow passing a heterogenous magnetic field has an anal- 
ogy with the flow around an obstacle. Increasing fur- 
ther the x-value of the vertical section and thus leaving 
the region with noticeable magnetic field, the Hartmann 
profile begins to transform back into a parabolic profile. 
Consequently, there must be a vertical movement from 
the bottom and top walls towards the center, see Fig. [8^ 
(note: we are only speaking about the vertical velocity 
component; the span-wise velocity component in section 
e is directed away from the center, compare with Fig. [Til. 
This phenomenon is opposite to that observed in Fig.[8]2: 
now there is no more braking force keeping the veloc- 
ity profile flattened in the vertical direction. As a result, 
the flow develops a swirling component in the channel 
corners. The swirling behavior extends far behind the 
magnet, see Fig. [8]/'. To our knowledge, this effect of 
swirling flow in the corners inside the rectangular chan- 
nel after passing the region of heterogeneous magnetic 
field has not been mentioned before elsewhere. 

To get more insight how the vertical velocity redistribu- 
tion takes place. Fig. |9] presents a vertical section along 
the midplane in stream-wise direction of the channel. 
The velocity vectors shown are the projection (m.v,Mz) 
of the total velocity. One clearly sees how the vertical 
velocity profile flattens before the magnet and becomes 
more convex after the magnet. Interestingly, the flatten- 
ing starts in the center of the profile because of the action 
of the Lorentz force while the decay of the Hartmann lay- 
ers begins at the top and the bottom wall because of the 
wall friction and viscosity in the fluid. The isolines of 
the normalized vertical profile p{z',x) in Fig. |9] give an 
additional information to judge how flat the profile is at 
a given stream-wise coordinate x. We define the normal- 
ized vertical profile as p{z',x) = Ux{x,0,z)/uo{x), where 
uo{x) = jf\ Ux{x,0,z)dz is the mean velocity in the ver- 
tical middle plane at given x and y = 0. Horizontal iso- 
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Figure 8: Vector plots of the vertical and span- wise velocity components at the following vertical sections of the 
channel: {a) x = —4, {b) x = —2, (c) x = 0, {d) x = 2, (e) x = 6, {f) x = \\ &t Re = 400, N = 36. An interesting 
feature is the swirling flow that arises in the comers behind the magnet (plot e,f). 



lines in Fig. |9] correspond to a profile, where the velocity 
changes in the vertical direction as for the parabolic pro- 
file in the inflow region, while the almost vertical isolines 
under the magnet represent regions with a flat profile, i.e. 
the Hartmann profile. Before (behind) the magnet the 
isolines converge (diverge) which demonstrates again the 
flattening (de-flattening) of the vertical velocity profile. 
The flattening process turns out to take place much faster 
(in four length units) than the de-flattening which is far 
from being finished in the six shown length units after the 
magnetic gap. 

4.2 Electric potential 

One main focus of this work is the demonstration that the 
flow under the magnet is determined by the value of the 



interaction parameter N, independently of the individual 
values of the Hartmann and Reynolds numbers as long 
as the Reynolds number is high enough. For this aim we 
first start with a comparison of two numerical computa- 
tions for the same interaction number N but with different 
Reynolds numbers. Fig. [TOl shows such a comparison at 
N =\6ior:Re = 900 (solid lines) and Re = 225 (dashed 
lines). One can actually see that the two contour line 
sets are not prominently distinct despite the fact that for 
both cases the Re numbers differ by a factor of 4, and the 
Ha numbers differ by a factor of 2. Before the magnetic 
gap the dashed and solid lines closely coincide with each 
other, and after the magnets there is a slight discrepancy 
since the inertial force for Re = 900 is higher than for 
Re = 225. More discrepancies one can find in the thick- 
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Figure 10: Electric potential distribution, in cen- 

tral horizontal plane (z = 0) for N = simulation: 
Re = 900 (solid) and 7? = 225 (dashed). 



Figure 9: Vector plot of the velocity components and 
u, in the central vertical plane at y = 0, /?e = 400, A'^ = 36. 
The isolines are plotted at constant level of the normal- 
ized vertical profile p{z\x) and give additional informa- 
tion about the flatness of the velocity profile. 



nesses of viscous boundary layers, however they are es- 
sentially of no importance since we are interested in the 
processes taking place in the central part of the magnetic 
gap- 
Now let us explain qualitatively the electric field distribu- 
tion in the system under consideration. In the quasi-static 
approximation the current according to Ohm's law eq. ^ 
is the sum of two terms: the electric field induced by the 
motion of the conducting fluid volume inside the mag- 
netic field, u X B, and the electric field E = — V(|) derived 
from an electric potential (|) generated inside the fluid 
volume because of the solenoidality of the electric cur- 
rent (due to Ampere's law V x B = /^j) and the isolating 
boundary conditions at the walls forcing currents to close 
inside the fluid volume. Taking the divergence of (O one 
gets the Poisson equation dUl for which the right hand 
side V • (u X B) = B • w plays the role of the inhomogene- 
itjH- Comparing with the Poisson equation for usual elec- 
trodynamics, A(|)(r) = — p(r) / (£()£), one sees that the ex- 
pression p* = — SqsB • w can be considered formally as 
an electric charge density induced by the movement of 
the electrically conducting fluid in the magnetic field. In 

^The second term u ■ (V x B) vanishes because the fluid volume 
is outside of the external magnet and induced magnetic fields are 
neglected in the quasi-static approximation. 



the duct flow, the largest contributions to B • w are gener- 
ated with the predominant vertical component of the 
magnetic field together with the span-wise velocity gra- 
dient dux/dy of the stream-wise velocity component, i.e. 
p* ^ z^zB^idux/dy). Therefore, looking in stream-wise 
direction, one finds a negative charge density near the left 
wall and a positive charge density near the right wall of 
the channel which correspond to the outer flanks of the 
M shape profile under the magnet. The electric field cre- 
ated by this charge density inside of the channel points 
in span-wise direction parallel to the y axi^ (The y 
axis points from the right to the left wall of the chan- 
nel). The braking action of the Lorentz force leads to a 
strongly diminished flow in the central region under the 
magnet and in this way the M shape profile is created. 
The inner flanks of the M shape profile provide opposite 
velocity gradients and opposite charge densities causing 
an electric field antiparallel to the j-axis in the centi^al 
region. It depends now on the strength of the magnetic 
field and the corresponding braking effect, whether there 
exists a region in the middle of the channel under the 
magnet, where the total electric field — the sum of both 
described charge densities from the outer and inner flank 
of the M shape profile — is antiparallel. In this case the 
electric field has to change its sign two times along the 
span-wise direction, corresponding to three sign changes 
for the electric potential. In the contour plot of the elec- 
tric potential, the stagnant region is manifested by closed 

"'Two remarks: (1) The electric field direction coincides, of course, 
with the direction of u x B. (2) Directly at the walls the charge 
density is high, but the induced electric field has to be zero because 
of the isolating boundary condition. 
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Figure 1 1 : Simulated span- wise electric potential profiles 
in the center of the magnetic gap (x = and z. = 0), Re = 
400 and = 4 (dot-dashed), = 9 (dashed), = 20.25 
(dotted), N = 36, experimental magnetic field (solid line 
1) and N = 36, simplified field (solid line 2). 



contour lines as we shall see below in Fig. [T3l [14] 
To illustrate the previous deduction about electric poten- 
tial behavior we plotted in Fig. [TT] span- wise profiles of 
(|) for different interaction parameter N. One can see that 
for the low interaction parameter N = 4 (dot-dashed) (^{y) 
shows clearly monotonic behavior, while for the high 
interaction parameter N = 36 (solid line 1) the sign of 
the electric potential changes three times. Correspond- 
ing curves in Fig.|6]and[TT]are plotted with the same line 
types. One clearly observes the tendency: the develop- 
ment of the inflection of the (^{y) curves (Fig. fTTI) is ac- 
companied by a lower minimum of the velocity profile 
(Fig. O around y = 0. 

The next two figures show how the electric potential 
changes by passing the magnetic field region. They 
are given in comparison with the experimental results. 
Fig. 1 12(a)] shows the profiles for N ^20 and Fig. |12(b)| 
shows the profiles for N ^9. As it should be, kinks of the 
electric potential in Fig ]12(a)] are more expressive than in 
Fig ] 12(b)] due to the larger interaction parameter. Also, 
the figures demonstrate that this effect is most expres- 
sive in the center of the magnetic gap {x = 0) where the 
magnetic field and, therefore, the braking Lorentz force 
is maximal. In the region of the inward (x = —2) and 
outward (x = 2) magnetic field gradient, electric poten- 



tial behavior is monotonic. The comparison between the 
simulated and experimental results reveals almost perfect 
accordance. 

To complete the comparison wit h the experimen- 
tal BAndreev. Kolesnikov. and ThessI (]2006f) 1 data we 
present also the contour plots. Fig. [13] shows the experi- 
mental and numerical results of the electric potential dis- 
tribution at N ^ 20 . Two contour plots give the level 
lines of the electric potential, the first, as it was directly 
and systematically measured in the experiment, and the 
second, as computed numerically. Qualitatively, the elec- 
tric potential distribution has the following main features: 
(i) it is an antisymmetric Al contour map with respect to 
the axis y = (ii) there are two global extrema directly 
at the side walls exactly under the peak of the magnetic 
field (x = 0) and (iii) there are two local sign alternating 
extrema slightly shifted in stream-wise direction at one 
third distance from the side wall (measured by the total 
span-wise width of the channel). Fig. [13] As one can 
see, all the features of the experimental electric potential 
measured at Re = 2000 are excellently reflected in the 
simulation at Re = 400 since both have the similar inter- 
action parameter, N 20. Another comparison for A/^ 9 
is given in Fig. [14] Now, the closed lines of the electric 
potential are weaker because of the lower interaction pa- 
rameter. The good accordance between experimental and 
simulated results is observed in Fig. [14] as well. 

Summing up, we find a very good agreement between the 
experimental and the numerical results comparing data 
for the same interaction parameter A'^. This holds even 
in complementary regimes with respect to the Reynolds 
number: turbulent inflow in the experiment, and laminar 
calculations in the simulation. This shows that the inter- 
action parameter indeed is the governing parameter for 
the flow under and near to the magnets, and that the flow 
in this region is determined by the magnetic and inertia 
forces. If the Reynolds number is not high enough, the 
viscous force begins to play a role as a third force in the 
system as can be observed to a small extent in our first 
comparison, see Fig.[lO] of this subsection. 

As the electric potential does not show a visible reaction 
on turbulence or no turbulence in the inflow, one could 
conclude that it is a quantity which is rather insensitive 
to different influences. We will show now that this is 
not the case, and consider for this aim the action of the 
span-wise decrease of the experimental magnetic field on 
the electric potential distribution. This dependence is 
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(a) (b) 
Figure 12: Experimental (symbols) and simulated (lines) span-wise electric potential profiles at z = and x = —2 
(crosses and solid lines), (dots and dashed lines), 2 (circles and dot-dashed lines). Plot a is for N ^ 20: {Re,N) = 
(2000, 18.6) (experiment) and {Re,N) = (400,20.25) (simulation). Plot b is for N^9: {Re,N) = (4000,9.3) (ex- 
periment) and {Re,N) = (400,9) (simulation). 



especially interesting, since many people believe that in 
most cases only the transverse component of the mag- 
netic field, changing along stream-wise coordinate, is of 
importance. For instance, most of the numerical simula- 
tions were performed only with B^(x) dependence, see, 
e.g. Jsterll (|l990h : Ulboussierel (|2004h l. To clarify the 
role of the magnetic field and demonstrate the sensitiv- 
ity of the <^{x,y) map, we fixed {Re,N) = (400,36) and 
compare two variants with the real experimental (show- 
ing j'-dependence) and simplified (span-wise uniform) 
magnetic field. As it is supposed typically, the simplified 
magnetic field has only stream-wise dependence in its in- 
tensity, and is constant at every fixed x in span-wise di- 
rection. The simulated results for the simplified field are 
presented in Fig. [15] One discovers that now the picture 
changes even qualitatively: there are no more closed lines 
of the electric potential, that is, the potential drops mono- 
tonically from one side wall to another (compare also the 
solid curves 1 and 2 in Fig. fTTl) . This is a strong indica- 
tion that the factors, which are decisive for the compar- 
ison of experimental and simulated results, are the max- 
imal interaction parameter in the magnetic gap and the 
proper configuration of the magnetic field. 



4.3 Electric currents paths 

Fig. [16] shows the closure of the electric currents in the 
central horizontal plane (z = 0) which were calculated as 
lines tangent at every point to the vector of electric cur- 
rent j. For sake of simplicity only few paths are shown, 
one in the region of increasing magnetic field, and others 
after the magnet in the region of decreasing field. In the 
case of constant magnetic field the loops of the electric 
current are located entirely in vertical planes and most of 
the current is concentrated in the Hartmann layers. By 
contrast, in the case of the heterogenous magnetic field, 
since an electric current intends to close itself in the re- 
gion where magnetic field is minimal (in order to close 
along paths with smallest resistance), one sees the cur- 
rent loops close themselves in the horizontal planes. The 
characteristic length of the loops is rather large as can 
be seen in Fig. [16] the turning point of j, starting under 
the peak of the magnetic field (x 0), is at \x\ f» 10. ..15, 
while a remarkable magnetic field intensity is felt only up 
to X 5. This fact must be taken into consideration when 
one selects where to put inlet and outlet of the numerical 
simulation box, otherwise it is possible to get artificial 
findings. 

Fig. [T6]reveals small closed loops for the electric current 
atx^ 2.2 and y ^ ±4.1. These loops are interesting be- 
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Figure 13: Electric potential distribution, ^{x,y), in the Figure 14: Electric potential distribution, (^{x,y) in cen- 

central horizontal plane (z = 0) for N ^20: experiment, tral horizontal plane (z = 0) for A/^ 9: experiment. 

Re = 2000, N=\S.6 (left) and simulation Re = 400, = Re = 4000,A^ = 9.3 (left) and simulation Re = 400, = 9 

20.25 (right). Bold dashed lines are physical borders of (right). Bold dashed hnes are physical borders of the 

the magnetic gap. magnetic gap. 



cause of the closure of the electric current is typically a 
big horizontal loop which envelops the space of the mag- 
netic field gradient. The geometric explanation for these 
small exotic loops is that they are projections to the XY 
plane from complicated electric paths developing essen- 
tially in 3D space. Fig. [l7] clearly shows these 3D paths: 
they are helices connecting the Hartmann layer and the 
middle horizontal plane. 

As well in our calculation presented in this paper, see 
Fig. [T7J as also in other calculations which are not shown 
here, helical current paths are only present when the level 
lines of the electric potential relief in the horizontal mid- 
plane contains closed lines, see Fig. [13] Our results 
strongly suggest that the presence of closed lines of the 
electric potential is a necessary condition for the appear- 
ance of helical current paths. As the span-wise uniform 
magnetic field does not lead to a triple change of the 
span-wise electric potential profile (i.e. no closed lines 
in Fig. [TSl) . the 3D paths of the electric currents for this 
case consequently are simple closures in the horizontal 
plane as is shown in Fig. [TSl 

There are also correlations between the span-wise inho- 
mogeneity of the magnetic field, 3D helical currents, and 
features of the stagnant region in M-shaped velocity pro- 
files. One sees a broad constant plateau in Ux{y) for the 



simplified field while for the experimental field there is 
a degraded minimum approaching zero (compare solid 
curves 1 and 2 in Fig. O. Since the level of the stagnant 
plateau for the span-wise uniform magnetic field is suffi- 
ciently high to keep electromotive contribution u x B in 
the electric current j, there is no double sign change of 
the electric potential (compare solid curves 1 and 2 
in Fig.fTTI). The stagnant plateau is not so broad but con- 
siderably lower for the experimental magnetic field, be- 
cause of the lower magnetic field strength near the side 
walls. The latter induces the braking Lorentz force near 
side walls to be smaller than in the center. As a result, 
the case of the real field looks somehow as a flow around 
the magnetic obstacle, while the simplified field is looks 
more similar to an uniform semi-penetrable barrier. 

The 3D helices of the electric current in the case of the 
real field (Fig. [TT] ) arise near the turning points where 
the electric field E changes its sign, hence £" ~ 0. These 
helices are maintained mainly by the electromotive com- 
ponent (u X B) in Ohm's law, eq. [3] On the other hand, 
these helices are located in the region of space where 
stream- wise component Bx{x) reaches the largest value, 
see Fig.m in the center of the outward magnetic field re- 
gion. Referring to the vector product u x B, we see that 
the term UyB^ is of importance for the size of the com- 
ponent, which is responsible for the uprise of the helix. 
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Figure 15: Electric potential distribution, ^ix^y) in cen- 
tral horizontal plane z = for /?e = 400, A'^ = 36, and 
magnetic field without span-wise dependence. 



Thus, in order to catch the helices one has to keep in the 
simulation all the components of the magnetic field. 
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Figure 16: Electric current paths in the central horizontal 
plane (z = 0) under the experimental magnetic field. Re = 
400, = 36. The lines are not equidistant since they 
were obtained by integration from a few starting points 
at {x,y) = (-0.5,0), (0.5,0), (2.1,4.1), (2.1, -4.1). 



5 Conclusions 



We have carried out numerical simulations for liquid 
metal channel flow under inhomogeneous magnetic field. 
For the computations the same channel geometry and 
the magnetic field configuration as in th e experiment of 
l Andreev. Kolesnikov. and ThessI ( 2006h l were used. 

Computations of the velocity field showed that all known 
general features — like for example the formation of 
Hartmann layers at the walls perpendicular to the main 
component of the magnetic field, or the formation of 
an M shaped profile in span-wise direction — are rep- 
resented correctly by our numerical code. All features 
of the flow together lead to a complex three-dimensional 
flow structure, which for the investigated regime of inter- 
action parameters, 4 < < 36, can only be determined 
numerically. As a new feature of the velocity field a 
swirling flow in the corners of the duct is observed. It 
begins shortly after the magnetic gap and extends far into 
the outflow region. 

The main goal of this work was to compare the electric 
potential distribution measured in the experiment with 
that of our numerical simulation. We were able to find 
a very good agreement of the electric potential distribu- 
tion for two sets of parameters which were used in the 
experiment ((i) = 18.6, Re = 2000, and (ii) A^ = 9.3 
Re = 4000) and the corresponding sets ((i) A^ = 20.25, 
Re = 400, and (ii) N = 9,Re = 400) in our numerical cal- 



culations. This comparison shows that the electric poten- 
tial distribution in the magnetic field region is solely de- 
termined by the value of the interaction parameter. This 
makes sense as in the magnetic field region the Lorentz 
force and the inertial force are strongest, while the vis- 
cous force is only important in the Hartmann layers and 
near to the side walls. The turbulence, which is present 
in the inflow of the experimental channel is in the mag- 
netic field region already negligible and has no influence 
on the electric potential distribution. A numerical test 
with a different inflow profile, which is more flat like a 
turbulent profile, also showed no difference. 

In addition, we could demonstrate that local extrema 
of the electric potential appear for interaction parameter 
higher than a critical value, which lies between 9 and 18. 
Simulations with a simplified magnetic field (only de- 
pendent on vertical and stream-wise coordinate, no span- 
wise field component) showed that local extrema of the 
electric potential map totally vanish in this case. This 
means that already a slight variation of the inhomoge- 
neous magnetic field can have a strong influence of the 
electric potential distribution. The same sensitivity on 
the magnetic field structure is found for the electric cur- 
rent density: In the case of the experimental magnetic 
field our computations revealed the appearance of com- 
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Figure 17: 3D electric current path for the experimental 
magnetic field. 



plicated helical current lines which are not present for the 
simplified magnetic field. 
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